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• Abstract 

^"1 \ The Wishart distribution on an homogeneous cone is a generalization of 

' ■^ | the Riesz distribution on a symmetric cone which corresponds to a given 

graph. The paper extends to this distribution, the famous Olkin and Rubin 
characterization of the ordinary Wishart distribution on symmetric matrices. 
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1^ : 1 Introduction 



, In many practical situations, there are manifest inter-relationships among several vari- 

ables. One important case is when several pair variables are conditionally independent, 
giving other remaining variables. For multivariate normal distribution, this corresponds 
to some zeros among the entries of the precision matrix. Due to this there has been an 
interest in distributions akin to the Wishart but defined more generally on various cones 
containing the cone of positive definite symmetric matrices as a special case. In partic- 
ular, Andersson and Wojnar [2] defined the Wishart distribution on homogeneous cones. 
This distribution is in fact an extension to homogeneous cone of the Riesz distribution on 
symmetric cones defined by Hassairi and Lajmi in [8]. In the present paper, we give a char- 
acterization of the Wishart distribution on homogeneous cones which is parallel to that 
given in [13] by Olkin and Rubin, or more generally in [6] by Casalis and Letac. In these 
papers, it is not assumed that densities exist, however distributions are assumed to be in- 
variant by the orthogonal group of the appropriate algebra. Our characterization uses the 
Laplace transform and a decomposition of a random variable on an homogeneous cone as 
a sum of random variables concentrated on certain subalgebras. The distribution of each 
component is then assumed to be invariant by the orthogonal group of the corresponding 
subalgebra. Let us give, a brief history of this characterization. Luckacs (see [11]) gave 
the following characterization of the Gamma distribution: Let X and Y be two non Dirac 
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and non negative independent random variables such that X + Y is positive almost surely, 
then U = X + Y is independent of V = X ^_ Y ^ an d only if there exist a > 0, p > 
and q > such that X and Y are distributed as Gamma distributions with parameters 
(p, a) and (q, a) respectively. Olkin and Rubin [13] extended this characterization to the 
Wishart distribution on the cone $7 of positive definite symmetric matrices. They showed 
that X and Y in £1 are Wishart if and only if (X + Y)- 1 / 2 X{X + Y)- 1 / 2 is independent of 
X + Y and its distribution is invariant by the orthogonal group K. The characterization 
has then been extended by Carter [5] to the Wishart on the cone of Hermitian matrices 
with entries in (D, and by Casalis and Letac [6] to any symmetric cone. There are other 
types of characterizations of the Wishart such as, for example, that given by Letac and 
Massam in [9], Geiger and Heckerman in [7], and Massam and Wesolowski in [12]. We 
also mention that more recently Bobecka and Wesolowski gave in [3] a characterization 
of the Wishart distribution on f2 without any assumption of invariance. However, they 
had to assume that the densities of both X and Y with respect to the Lebesgue measure 
exist and are twice differentiable. The extension of this characterization to homogeneous 
cones has been the object of a paper by Boutouria [4]. The paper is organized as follows. 
In §2, we recall various definitions and preliminary results relevant to the Wishart on 
homogeneous cones and we establish some results concerning the Vinberg multiplication 
and determinant calculation. In §3, we state and prove our main characterization result. 



2 The Wishart distribution on homogeneous cones 

For the convenience of the reader, we will give here the elements of Vinberg algebras 
and homogeneous cones essential to working with the family of Wishart distributions on 
homogeneous cones. These elements are taken from [2] and the reader is refereed to this 
paper for further details. After recalling this development, we give three examples. For 
the third example, the poset is isomorphic to a rooted tree while this is not so for the 
second example. The first example corresponds to the ordinary Wishart distribution on 
symmetric matrices. 

Let 7 be a partially ordered finite set (herewith abbreviated as poset) equipped with 
a relation denoted ■<. We will write i -< j if i < j and i ^ j. We assume that / satisfies 
the following condition 



(F): 



for any two points i and j in / such that either i < j or j < i 
the path on the Hasse diagram of / between i and j is unique. 



For all pairs £ I x I with j -< i, let be a finite-dimensional vector space over IR 
with riij = dim(Eij) > 0. Set 



Aij 



IR 

E, 



Eji 
{0} 



for 
for 
for 



i=3 
j < i 
i < j 
otherwise 



and A = \\ A 

ijeixi 

iei 



n ■ 



n 



Define n,- 



n.i 



1 + -(raj. + n.i), i G / and 
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An element A = (a^, i,j G I) of A may be seen as a matrix and so we define its trace as 

tiA = y~]au. 

iei 

Let fij : Eij — > Eij, i y j, be involutional linear mappings, i.e., f^ 1 = fij. They induce 
an involutional mapping ( A i->- A* ) of A given as follows: A* = (a*j\(i,j) G I x I), where 

{an for i = j 

fij(aji) for j ~< i 

fji(aji) for i<j 

otherwise. 

We now define the following subspaces of A: 
the upper triangular matrices 

T u = {A = ( aij ) €A,Vi,j€l:i£j=> Oij = 0}; (2.1) 
the lower triangular matrices 

Ti = {A = fa) eA,Vi,j el :j j<i^> a tj = 0}; (2.2) 
and the Hermitian matrices % = {A £ A, A* = A}. 

The sets of upper and lower triangular matrices in A with positive diagonal elements are 
respectively denoted by and 7^ + . The sets of upper and lower triangular matrices 
with all diagonal elements equal to 1 are respectively denoted by and 7^. The sets 
of diagonal matrices and of diagonal matrices with positive entries are denoted by V and 
D + , respectively. 

We are going to equip the vector space A with a bilinear map called multiplication and 
denoted by (A,B) >->■ AB. For this purpose we need to define bilinear mappings Aij x 
Ajk — > Aik, denoted by (aij,bjk) dijbjk-, and then define AB = C = (cij\(i,j) G I x I) 

ixei 

The multiplication is required to satisfy the following properties: 

i) \/A G A; A / tr(AA*) > 
it) VA,B G A; (AB)* = B*A* 

Hi) VA,B G A; tr(AB) = tr(BA) (2.3) 

iv) MA, B,C eA; tr(A(BC)) = ti((AB)C) 

v) VC/, S,TeTf, (ST)U = S(TU) 

vi) VJ7, T G 71; T(UU*) = (TU)U*. 

An algebra A with the above structure and properties is called a Vinberg algebra (For 
more details, we can refer to [2]). We choose the element A = (aij\(i,j) G I x I) of V such 
that an = 1, V i G I as the unit element of „4 and we denote it by e. Vinberg proved in 
[14] that the subset V = {TT* G A, T £ T^} CH C A forms a homogeneous cone, that 
is the action of its automorphism group is transitive. 

The definition of V could be changed to the following equivalent definition 

V = {TUT € A, T £ TJ 1 , D G £>+}. 



3 



The two decompositions S = TT*, T G % + and 5 = T x DTf, T\ G Tj 1 , D G V + are 
unique and their connection is given by T = T\^f~D where = diag(\/c#i, i G /) G f + 
when D = diag(<ij, i G I) G £> + . 
For S 1 = (sy, i,j e I) = T X DT*, we write Ai =<%].• 

If ^ opp is the opposite ordering on the index set I, i.e., i ^ opp j j ^ i. The 
Vinberg algebra A opp = -4i°; PP > where 

i,j£lxl 

IR for i = j 

for j -< opp i 
£y for i -< opp j 
{0} otherwise, 

differs from the Vinberg algebra A only in the ordering of the index set I. It is proved 
(see [12]) that o PP = {T*T G A, T G 7^ + } is the dual cone of V which is also denoted 
V* . 

Let's explain why condition (F) on the poset I is required for the definition of a Vinberg 
algebra. In fact, the property vi) in (|2.3p fails to be verified if the condition (F) is not 
satisfied. Suppose that (F) is not satisfied, then there exist two elements i and j in / such 
that either % ~< j or j ~< % and the path on the Hasse diagram of / between i and j is not 
unique. Without loss of generality, we suppose that i -< j. Then there exist k and s in 
/ such that k^s, i~<k~<j and i ~< s ~< j. Consider the elements T = (t nm ) njme i and 
U = (M nm ) n , rae j of 71 defined by t nm ^ and u nm ^ if n,m G {i,j, k} and t nm = and 
Unm = otherwise. Then 

[T(UU*)]jk = tjiUuUik + tjk(uf k + u 2 kk ) + tjj{ukiUij + UjkUkk) 

and 

[(Tf7)£/*]jfc = (tjiUu + tjkUki + tjs^si + tjjUji)Uik + (tjfcUfcfe + tjjUjk)Ukk- 

We see that T(UU*) ^ {TU)U* so that m) in J33D is not satisfied. 

Andersson, Letac and Massam [1] have also shown, without using the theory of Vinberg 
algebras, that if (F) is satisfied then, the cone V is homogeneous. 

Let G be the connected component of the identity in Aut("P), the group of linear 
transformations leaving V invariant. We recall that % '■ G l— * ^+ i s sa fo to be a multiplier 
on the group G if it is continuous, x( e ) = 1 an d xihihi) = x(^-i)x(^2) for all hi, hi G G. 
Consider the map vr : Te7j + 4 vr(T) G ^(T^) C G such that for X = WW* G V, 
WET+ 

tt(T)(X) = (7W)(W*T*). (2.4) 

Andersson and Wojnar show in [2] that the restriction of a multiplier x to the (lower) 
triangular group 7j , i.e., pir : TJ + — > M+ is in one to one correspondence with the set 
of G I) G 1R 7 . To each % corresponds a unique (up to a multiplicative constant) 

equivariant measure v x on C with multiplier x under the action of G, that is for all h G G, 
the image measure h~ x v x of z^ x by h^ 1 is 

h~ l v x = x{h)i> x . 



A opp 
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For G V* , the Laplace transform of u x is 

L x {9) = I exp{-9(P)}du x (P). (2.5) 
Jv 



We define the ^-inverse of 9 by 



9 x = -^\ogL x {9). 



The mapping 

is a bijection. Its inverse is denoted by 

FeP i — y P~ x eV*, 

and we have the properties 

(h^P)-* = l hP- x and ( t h9) x = h~ 1 (6 x ), (2.6) 

where h £ G, P £ V and 6 £ V* . Also for G "P*, it is convenient to introduce the 
notation a = 9~ x . Then (|2.5p may be written as a function of a as 



n x (a) 



■ [ exp{-a- x (P)}dv x (P), 
JV 

so that for h G G, a G "P, n x has the property 

n*(/w7) = X (h)n x (a). (2.7) 
Andersson and Wojnar (see [2]) consider the set 



X = { x o 7T : 7^+ -> R+; Ai > -y , i G J}, 



and show that for x £ ^ the measure concentrated on V 

u x {dX) = l[x^ n n v (X)dX (21 



iei 

generates the Wishart natural exponential family of distributions on V absolutely contin- 
uous with respect to the Lebesgue measure, parameterized by a G V 

^ n ^ n *5p 

HW x , a (dX) = ; eJ n / eJ Xi exp-{ tr (a^X)}l v (X)dX. (2.9) 

11 1 (A* - -a 1 -) 11 ^m. 

It is shown in [2] that if a = ZZ*, where Z G 7^ + , 

cj-x = (Z*)" 1 diag(A i |i G t)Z~ x . (2.10) 

We are now going to give three examples of homogeneous cones and their corresponding 
Wishart distributions. 
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Example 2.1 In this example, we show how the ordinary Wishart distribution on sym- 
metric matrices is a particular case of the general Wishart distribution HW X)CT on homo- 
geneous cones. Take I = {1, ••-,/} equipped with the usual total ordering X on integers, 
denote by I its cardinality, and set Eij = M, j -< i. Then the vector space A is the 
space M(I, M) of all I x I matrices. With the standard multiplication and inner product 
the vector space A4(I,M) is a Vinberg algebra. The homogeneous cone V in this Vinberg 
algebra is the usual cone of I x I positive definite symmetric matrices usually denoted 0. 
In this case, we have 

n.i = I - i, ni.=i- 1, ni = — — , re.. = 

and every multiplier \ : G — > has the form x(h) = \ det(h)\ a , a £ M. Thus for all 
i & I, Aj doesn't depend on i, it is equal to A = ^Qot, so that x can be replaced by A. 
From \2.10\) . it follows that a~ x = Act -1 , where cr _1 is the usual inverse of a in V . Since 
for S = (sij, i,j £ J) £ fl, det(iS') = JJs^]., the density of the Wishart distribution \2. 9\) 

iei 

becomes 



W x> a (dS) = — — exvi-Xtria-^yiaWdS, (2.11) 

iei 

which is the usual Wishart distribution on Q with parameters A > and a & tt. 
Example 2.2 Let I = {1,2,3,4} be an index set and consider the poset defined by 

1 -< 3,2 -< 3,2 -< 4. 

The homogeneous cone corresponding to this poset cannot be of the type Qq of incomplete 
matrices with submatrices corresponding to the cliques of G being positive definite. We 
will not give the details of the argument for this assertion here. However, an accurate, 
albeit short, argument is that according to Theorem 2.2 in [10], this can only happen if the 
undirected graph obtained by dropping the directions on the graphical representation of the 
poset is a homogeneous graph. The graphical representation of the poset in this example 
is the directed graph with directed edges 

{(1,3), (2,3), (2,4)}. 

The undirected graph obtained by dropping the directions is the three-link chain with undi- 
rected edges {(1, 3), (3, 2), (2, 4)}, which, by definition, is not a homogeneous graph. 

Let us illustrate this in the case where = IR for all appropriate (i,j). The Vinberg 
algebra A is thus the vector space of all I x I real matrices with zeros at the (1,2), (2,1), 
(3,4) and (4,3) entries. An element X = (xij)ij e j of V is decomposed into X = TDT* 
with T € 7j j D G T>, such that 

ld x \ 
d 2 
d 3 ' 
V d 4 J 



/TJ* 



/l (is \ 

1 *23 *24 

1 

Vo II 



,D 
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Then the equation 



X 



1 Xi 





XYZ 


^ 







X2 


X23 


2?24 




X31 


X32 


X3 







\ 


X42 





x 4 J 




/ dl 












0I2 




^2*23 



*31^1 ^32^2 *32 0^2*23 + ^3 
^42^2 t 32 d2t23 + d 3 



\ t 42 d 2 

/ias a unique solution 



£42^2*23 



/ 1 0\ 

10 

*3i *32 1 

V t 42 1 J 
\ 

^2^24 
^32^2*24 
^32^2*24 
^42^2^24 + ^4 / 



/di \ 
d 2 




\0 



d 3 
d 4 / 



/l iis \ 

1 *23 *24 

1 

\0 1 / 



di= Xi= ii 3 



3 = x 3 - x 32 x 2 X 23 



g!3 
X\ 

X[ 3 ]., 



d 2 = X 2 = X[ 2 ]., i 2 3 = 
d 4 = X 4 — X 4 2X 2 ~ 1 X24 



2?23 



X2 
= X [4 ].. 



, *24 



X24 
^2 : 



(2.12) 



TTte cone is then the set of symmetric matrices X = (x^) swc/t that 
x\ > 0, X2 > 0, X3 — x 2 ~ 1 x 2 3 > and x 4 — x 2 " 1 X2 4 > 0. 
Since Ey = M, we have that riij = 1, n±. = 0, ra 2 . = 0, 723. = 2 and n 4 . = 1. Therefore 

X = {(Ai, A 2 , A 3 , A 4 ) : Ai > 0, A 2 > 0, A 3 > 1, A 4 > -}. 
and the Wishart distribution on this homogeneous cone is given for and a G V by 

A Q 4 



^ n A^r f x*rVjr 2 xH 



W X)tr (dX) 



1=1 



'[1]. ^[2]. -"[3]. -"[4]. 



r(A 1 )r(A 2 )r(A 3 - i)r(A 4 - ±) n 



■exp-{ tr (CT" x X)}ip(X)dX 



i=i 



Example 2.3 In i/izs example, the poset is isomorphic to a rooted tree and therefore the 
homogeneous cone V is of the Qg type, that is, it corresponds to a homogeneous graph 
(see [10], Theorem 2.2). 

Let the index set be I = {1, . . . , k} and the poset be defined by 

1 -< i, i = 2, . . . , k. 

We can immediately see here that the undirected graph obtained by dropping the directions 
on the directed graphical representation of the poset is the star-shaped graph with k vertices 
and undirected edges i = 2, . . . , k}, which is a homogeneous graph. According to 

Theorem 2.2 in [10], the cone V is therefore of the Qg type. Again, we take Eij = M for 
all appropriate (i,j). Then n_± = k — 1, n\. = 0, n± = and n.i = 0, n^. = 1, ni = |, 
Vi / 1. The Vinberg algebra A is the vector space of all A = (a^ £ / x I) G M(I, M) 
with a,ij = when i and j are not related. We have 



X = {(Aj, i = l,...,k): Ai > 0, A* > 



k}, 



and the Wishart distribution is given by 



ii aNJj. 2 ii 4* J 



HW x , a (dX) 



i=l 



i=2 



r(Ai) n r(A, - 1) n < 

i=2 1 



exp-{ tr 



A, 



i=l 



where x = {A^, z G 1} G <Y and a x is defined as in \2.10\) . for a G V ■ 

To illustrate the decomposition of an element of V in this example, we consider the case 

/ Xl Xu Xi3 X U \ 

where k = 4. An element X = X21 X ^ ^ _ 

x 31 x 3 

\ X41 X4 ) 

incomplete matrix 



of V can then be written as an 





/ Xl 


X12 


X13 


xu\ 


X* = 


X21 


x 2 


* 


* 




x 3 i 


* 


X3 


* 




\XAl 


* 


* 


x 4 / 



It can also be decomposed as X = TDT* with T G 77 , -D G X>, where 



/I 


tl2 




tu\ 





1 














1 













1 / 



D 



( di 
d 2 







d 3 
\ d 4 J 



The equation 



X 



( 1 


0\ 




(d 1 





\ 




(1 ii5 


*13 


tl4 


\, 


t2l 1 










d 2 







1 










isi 


1 







d 3 










1 







W 


l) 




\o 





dj 




\0 





1 


J 


/ di 


ti2di 




t 13 d 


1 


ii4<i 


i ^ 






£21^1 


t2ld\t\2 


+ d 2 
















^31^1 









^31^1^13 


+ d 3 











\ti\dx 














t4id\tu 


+ d 4 J 







has a unique solution 
d\ = x\ 



X\2 , X 13 X U 

X\i]., I12 — ,tl3 — , tl4 — 

1 ' Xl X\ Xl 



do 



X2 - X2lX 1 X12 



X[2]- ,d 3 = X 3 - x 3 ix 1 1 xi 3 



X\ — X4lX 1 X\H 



Hence, in this example, the cone V is the set of symmetric matrices X such that 

xi > 0, X2 — x\ X x\ 2 > 0, x 3 — x^ x x\ 3 > and X4 — x^ l x\^ > 0. 

Let us note that the distribution of X £ Qq can be obtained from the distributions of 
di,i = 1, . . . , k and tu, i = 2, . . . , k. This transformation is given explicitly in [10], Formula 
(3.15) and Theorem 4-5. 
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Now, we recall that an edge is written as an arrow from its origin to its destination, 
and we define the in-degree of a vertex to be the number of edges having this vertex as 
their destination. A vertex is considered a source in a graph if its in-degree is (no vertices 
have a source as their destination) and if it is the parent of at least two vertices. This 
notion is needed in the following result. 

Proposition 2.1 The Vinberg multiplication TT* is the same as the standard matrix 
multiplication T.T* if and only if there is no source in I. 

Proof (=>) Suppose that there exists a source a in /. Let b and c be in / such that 
a -< b, a -< c, and consider T = (%)ije/ G T t + such that % c = 0. Then T.T* ^ TT*, 
which contradicts our assumption and therefore the source does not exist. 
(<=) Suppose that there exists T such that T.T* ^ TT* . Then there exist in I two 
different and non connected element i and j such that {TT*)ij = and (T.T*)ij ^ 0. As 



In what follows, we suppose that I has no source and we call m a maximal element in /, 
if it is not less than any element of /. 

Proposition 2.2 Let m be a maximal element in I and let J-< m = {i I E I, % ^ m}. 
Then has no source. 

Proof Suppose that there exists a source a £ /-< m and let b and c in be two children 
of a with a < b, a -< c. As m is a maximal, we have that a ^ m. Also m ^ b, because, if 
m = b, then c^b. Thus b ■< m, c ^ m and a ■< m. Hence we obtain a -< b, a -< c, b ~< m, 
c -< m and a -< m which is impossible by Condition (F) of this section, and the lemma is 
proved. □ 

For T = (tij)ijei is in T t + , we define the element T^ of Ti by T^ = (t^)j je /, with 
t'jk = tjk if i ^ h k and t'- k = otherwise. It follows from Proposition 12.21 that if m is 
maximal, then for any Z 6 V, we have 



Indeed, m maximal implies that I~^ m has no source, and in this case, the Vinberg multi- 
plication is nothing but the standard multiplication of matrices. 
Also, we have that for any i E I, 




& is a source which is a contradiction. 



□ 



det Z^yn 



n %' 



det Z^i 



jQ Zjj/j. and det Z <{ = Z^_. 




because i is maximal in 1^. Therefore, for all v £ I, 



det Z~( 



(2.13) 



bi- 



det Z^ 
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We now introduce the notion of a division algorithm in a homogeneous cone needed 
for our characterization of the Wishart distribution on the cone V . This notion is defined 
in [6] in the case of a symmetric cone. A division algorithm is a measurable map 

9--V -> G 

U i y g(U) 

such that g{U){U) = e. In particular, for U = TT* , we define g(U) = ir{T~ l ), where it is 
defined in ()2.4p . then g is a division algorithm, so that if X = WW* G P, then 

5 (*7)(X) = (T~ 1 W)(W*(T~ 1 )*). 

This algorithm is the one that we will use in all what follows. 



3 Main characterization result 

In this section, we state and prove our main characterization result concerning the Wishart 
distribution on homogeneous cones in the line of the characterizations given for the ordi- 
nary Wishart on symmetric matrices by Olkin and Rubin [13] and by Casalis and Letac 
[6]. Our considerations here will be restricted to the case of homogeneous cones with 
E ij =TR, (ij)elxl. 

It is easy to see that the Laplace transform of the Wishart distribution HW x<a on the 
homogeneous cone V is given for 9 G V* by 

n x[(g + a -x)x] 

L ™ X A0) = — — • 

As n x (o~) = J^°"^j n x (e), this can be written as 

Lhw x ,(t(6) = - eI 



iei 



II 



iei 



Recall that for i £ I, we denote 1^ = {j G I; j ■< i} and = {j G I; i ^ j}. 
Let i\ and 12 in /, such that i\ ^ 12 and 12 ^ ii, we say that j separates i\ and 12 if 
j G n Ij 2 -< and j {ii, 22}- In this case, j is called a separator. We also denote 
S{ = {j G Ij-<; J is a separator and V 7^ j, k -fc j}, 6 = Si and 5 = {i G ©, V j 7^ 

i, j -fi i}. This leads to the following decomposition of an element of V which will serve 
in our characterization result. 

Proposition 3.1 Let Z = TT* be an element of V with T G T t + ■ Denote p = {i G 
/, I^i = 0} and define, for i E I, 



T-<TL and Zi 



Zi< - Z s ^ if i G p 

seSi 

Zi^ if i G S 

otherwise, 
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Then we have that 



Z = E Z * = E z - 



(3.14) 



Proof: In order to prove the equality (I3.14p . we compare the blocks of Z and Z% 



on each subalgebra Aj = Am, j G I which we denote respectively by (Z)j and 

( Zi)j. From the definition of (Zj^), we first observe that 



(Z)j - {Zj^)j, 



V jel. 



(3.15) 



Now we discuss according to the position of j. 

If j G" p U S, then there exists a unique io G p U 5 such that j G ^i d- The fact that 
j G Ji -< implies that (Zi )j = (Zj^)j. From this and (|3. 15)) . we get 

(Z)jf = (Zj-<)j = {Zi )j = ( ^ 

iepus 

If j G S 1 , we have {Z{)j = 0, for all i G p. Therefore 



For j G p, we have that ( ^"] Zi)j = Zj = Zj^ — Z s ~<,, and we consider separately the 

iGp seSj 

cases Sj = and Sj ^ 0. If 5j = 0, then (V^ = and (Z)j = (Zj-<)j = Zj = ( ^""] -^j)j- 

ieS iGp 

If 5,- / 0, = (E = (E Z *^)i- Therefore (^)i = ( E Z i)r D 

iGS ie5j ie5j iGpUS 

Example 3.1 Consider the following poset on I = {1,2,3,4} where 1 -< 3, 1 -< 4, 2^4. 
T/zen 5 = Si = 5 2 = {4}, J 3 ^ = {3}. Hence Z 3 = (0), 



( \ 





\ t u j 



( \ 





V £ 4 4 J 



and 



Z, 



( *n \ 



tl3 i 3 3 

V ti4 o o t 44 / 



/ \ 

t 22 o 



V t 24 t 4 4 / 



/ in t 13 t u \ 



o o t 33 o 
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Consider the Vinberg subalgebra of A defined by Ai = Y\ ^-kh an d denote Vi, Gi and 

ei respectively, the corresponding homogeneous cone, connected component of the identity 
in Aut(Vi), and unit element. Also denote K = {g £ G, g(e) = e} the orthogonal group 
of A, and Ki = {k £ K, k(e{) = e{\. Finally let 

gi. Vi -> d (3.16) 

such that 9i (U^){X idl ) = (TrjW^)(W^(Trj)*), where = W^W*^. Then it is easy 
to see that 

(g(X + Y)X)^ = 9i {X^ + Y^)X^- 

From now on, a Wishart distribution HW x ^ e x will be called a standard Wishart dis- 
tribution. Next, we verify that any Wishart distribution on the homogeneous cone V may 
be standardized by a linear transformation. 

Proposition 3.2 Let X be a random variable valued in V. Then X is HW Xt(7 if and only 
if there exist p in G such that p(X) is HW Xte x ■ 

Proof (=>) Suppose that X is HW X ^ and write a = TT* with T £ , x = {Ai, i € 
I}. We have that y / diag(Aj, i £ I)T~ l € 77 , we then consider the element of G, p = 
7r(Y / diag(Aj, i G I)T~ l ), where tt is defined by (|2.4h . Using (|2,6p and (|2.7|) . we have for 
9 e V*, 

L P{X) (9) = E{e-^ X ))) 

Ul(p*(9) + o-x)^l 
iei 



iei 



\{[p-\e+p*-\a-^r]^ 



iei 



n 

iei 

and using the fact that n x (g^ 1 (u)(a)) = x{d~ l {u))n x (a) and x{d{ u )) = rj u [i] Al ( see [2])> 
we easily verify that 

n^>)^i.=n^i. 

i&I iei 

This applied to = p gives 



l p(x)(Q) 



iei 

in*** 



iei 



12 



As p* 1 (a x ) = )) x = (e x ) x = e, we obtain 



iei 

which is the Laplace transform of a HW Xte x distribution. 

(<=) Suppose that there exist p in G such that X' = p(X) is HW x ^ e x. Then using again 
(USD and ||2ZD, we have for 6 G p*, 

L x (0) = E(e- tT ^) 

HKp*-\e) + e)x}*i 

iei 

iei 

iei 



Hence X is tfW^, with er = (p*(e))* eP. □ 

We are now in a position to give our characterization results. According to Proposition 
13.21 the statements will concern the Wishart distribution HW x ^ e x. 

Theorem 3.1 i) Let X = TT* , with T = (^)i,ieJ *n ', be a random variable with 
Wishart distribution HW Xy(J . Then o~ = e x if and only if the tij, i, j G I are independent, 
ii) If X and Y are two independent random variables with respective Wishart distribution, 
HW X:e x and HW X , eX > , then for i G I, the distribution of = <?j(JQ^ + Y^X^ is iQ 
invariant, where gi is defined by \3.16\) . 

Next, we give the reciprocal of this theorem. 



Theorem 3.2 Let X and Y be independent random variables valued in V . Write X = 
TT* and Y = MM* , with T = (%)ijeJ and M = (my)ij- e j in T t + . Consider the divisions 
algorithms g defined by {2.1J$ and gi defined by 13.16]) and suppose that 

(i) the tij, i, j G I are independent and the rriij, i, j G I are independent, 

(ii) X + Y is independent of g(X + Y)(X), 

(Hi) for i G p U S, the distribution of = gi{X^ + Y^)Aj^ is Ki invariant. 
Then there exist x, x' £ % such that X ~ HW x ^ e x and Y ~ HW X , eX i . 

Before embarking on the proofs of these theorems, it is worth mentioning that in the 
particular case where the Vinberg algebra is the algebra A4(I, H) of I x I matrices (see 
Example 12. ip , Theorem 13.11 and Theorem 13.21 give together the famous Olkin and Rubin 
characterization of the ordinary Wishart distribution on symmetric matrices. In fact 
in this case we have, p = {!}, S' = 0, so that p U S = {!}, and it follows that only 
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K\ = K = {g £ G, g(e) = e} appears in Theorem the point (iii) of 13.21 We also have 
e x = Ae. Hence Theorem 13 . 1 1 b ecomes : 

i) Let X = TT*, with T = (iij)jjeJ m ■> De a random variable with Wishart distribution 
W\,(t- Then a = Ae if and only if the tij, i, j E I are independent 

ii) If X and Y are two independent random variables with respective Wishart distribution, 
W\ t x e and W\' t \i e , then for i £ I, the distribution of gi(X^ + Yi^jX^ is K{ invariant, 
where gi is defined by (|3.16p . in particular g(X + Y)(X) is K invariant. 

For the proof of the theorems, we need to establish the following result. 

Lemma 3.1 Let X and Y be two independent random variables with respective Wishart 
distribution, HW Xte x and HW X , eX i . For i £ I, let ctjj.^ denote the distribution of X^ 
conditional on X^ + Y^ = U^. Then for f in Gi, the image measure fctxj^ of ajj^ by 
f is such that fotjj^ = OL^jj.^y 

Proof Let H : A. — > H and F : A — > 1R be any continuous functions with compact 
support. Let Xi an d x'% De t wo multipliers on Gi and set xj = XiXi- Denote v Xi , u x < and 
v x i the corresponding equivariant measures concentrated on Vi as defined in (I2.8f) and set 
X\^ = f(Xi^) and Yjz, = fiY^). We consider the following equalities obtained by using 
the decomposition 

^(dXi^idY^) = au^dYi^idU^) 
and using the change of variable formula for integrals. For = X^ + Y^, we have 

J = [ H(Uii)F(XlJf( aUi J(dXlJi*\dUii) 
J A 2 ~ 

= [ H(U^)(Fof)(X^) au ^(dXi^) V x '{dU ir ,) 
J A 2 

= [ H(Xi^+Yi^)Fof(Xi^(dXi^(dY^y 
J A 2 

Thus 

J = [ Hif-HxU + Y^FiXUXfu^idXUXfu^dY^ 
J A 2 - - - 

= Xiir 1 ) I H{r\xU+Y^))F{xU)u^{dxU)v<{dY> < ) 
J A 2 - - - - - 

= Mr 1 ) f Hif-HS^mx^as^dX^idSi^) , 
J A 2 

where = Xf^, + Y-z,. Since f~ 1 (Si^) = U^, we have 

J = xHr 1 ) [ H{U^)F{XU)a m ^{dXU)r^ x hdU^) 
J A 2 

= [ H(U^)F{xU)a m<) {dXU)v x '{dU^y 
J A 2 

Comparing this expression with the definition of J, the lemma is proved. □ 

Proof of Theorem 13.11 i) (<=) See [2] where it is proved that if X ~ W X;e x, then 
%> hj G I are independent. 
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(=>) Suppose that X = TT* , with T = (ty)<je/ in 7 i + ' is HW x,° with a eX - We 
will show that in this case, the ijj are not all independent. As a 7^ e x , there exist 
i j such that erjj 7^ 0. Writing cr = WW*, with W = (ufyOtjeJ S 7^ + , we have 
that (Tjj = WikWjk, and as cr^ 7^ 0, there exists k £ I such that Wi k Wj k 7^ 0. Prom 

fee/ 

Proposition 13.21 we have that X' = vr( v / diag(A i , i G J)VK _1 )(X) is #W Xie x. If we set 



X' = 55*, with 5 = (sij)ij^i £ T t + and Z = TPydiag(A i 1 , i £ I) = (zij)ij e j, we can 
write 



X = Tr-Hv/diagCAi, i i /JW- 1 )^') = (ZS)(S*Z*). 
This implies that tj^ = Zjp s iMk f° r J 7^ Hence 



fJ.€l 



E(tjktkk) — E(tjk)E(t kk ) — Zj fl z kk E(s ll kSkk) — ^ z jn z kkE{s flk )E(s kk ) 

Ate/ Ate/ 

= z kk z j^E(s^Skk) + z jk z kkE(s\ k ) — z j^ z kkE(s^k)E(Skk) ~ z jk z kkE(s k kY 



fij^k /j,^=k 



We use the fact that the Sjj are independent, because X' is HW x ^ e x, to obtain 
E{tjktkk) — E(tjk)E(tkk) = z jk z kk{E{s\ k ) — E(skk) 2 ) 

= Zj k Z kk Y&x{s kk )- 



This is different from because 2^ = w kk \ A., 1 7^ 0, z,-jj. = Wj k \ \ hk 7^ and s^fc is not 



degenerate, since s? fc is gamma distributed (see [2]). 

ii) Using the notation of Lemma 13. II and the fact that U and V are independent, for two 
arbitrary continuous functions with compact support, H : A — > H and F : A — > 1R, we 
have that 

E{H{U^)F{V^)) = E{H{U^) [ F{ 9i {U id )X ldi )a Ux< {dX id )). 

J A 

Writing = gi(Ui^,)Xi^ in the last integral and using Lemma 4.1, we get 

E(H{U^)F(V^)) = E(H(Ui±) [ F(V^)a ei (dV^)) 

J A 

This proves (ii). □ 

Proof of Theorem 13.21 Let i € I such that X{ 7^ 0, where is defined as in (|3.14p and 

let U = X + Y = WW*, V = g(X + Y)(X). Without loss of generality, we can suppose 
that 5j has just one element; 5j = {s}. We will consider first the case where i £ S. In this 
case, from the hypotheses of independence, we have, for A\, A2, B\, B2 and C\ in V*, 

B(exptr(AiWi d - A 2 W S ± + BxW^W^ - B 2 W S ±W^ + CiV^)) 

= E(exptT(AiW^ - A 2 W S ^ + BAV^WU - B 2 W s ^W^))E(exptv(C 1 V i ^)) 

= f i (A 1 ,A 2 ,B 1 ,B 2 )h i (C 1 ), (3.17) 
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where fi(A 1} A 2 , B 1} B 2 ) = £(exptr(,4iW^ - A 2 W a ± + B x W id WU ~ B 2 W S ±W S %)) 
and hi(Cx) = £(exptrCiV^). 
We adopt the notations: 

and we define, for B G V* , fi(B) = E(exptrB(X + Y)i)- 
It is clear that if A\ = A 2 = and B\ = B 2 = B, then 

(fi)jk = (fi)jk when {j,k}<£Si (3.18) 

and 

(fi)jk,in = {fi)jk,ln when {j, k} (£ Si and {l,n} <£ Si. (3.19) 

Since we can differentiate under the expectation, there is a relation between the second 
partial derivatives with respect to A\ and the first partial derivatives with respect to B\, 
namely we have 

V 92 - g dQ f> f32 o) 

^ d(A 1 ) ja d(A 1 ) ka dh--. dt q d(B 1 ) jk dh---dt q ' { ■ ' 

where t\,- ■ -,t q are any arguments of /j. 

Differentiating fj3. 171) successively with respect to (Ai)j\, (Ci)^, and (Ai)^, for 
{ j, k} <£. Si and summing over A and [A, give the basic differential equations 

E{{Xi) 3k exptr^W^ - A 2 W S ± + BtW^W*, - B 2 W sd WU + C%V^)) 

= E // A - A " j/ '').V- ( 3 - 21 ) 

And differentiating (|3.2ip successively with respect to (Ai)i u , {C\) va and (Ai) na , for {/, n} (£ 
Si and summing, we obtain 

E((Xi) jk (Xi) ln exptr(AiWi-< - A 2 W S ± + B X W^WU - B 2 W sd W^ + CM*)) 
= £ ft^^^x^ (3.22) 

Now, we use the hypothesis of invariance of the distribution of Vi-< by the orthogonal 
group Ki. In terms of Laplace transform, we have that for any k G Ki, 

hi(C x ) = £(exptrCi^) = £?(exptrCiA;(Vi-<)). (3.23) 

By the choice of suitable k in Ki and differentiation of (|3.23p . we first establish that there 
exists a real constant 0i such that 

dh (C ) 

(hi) jm (0) = * 1 | Cl=0 = 0iS jm , forj, meli-K, (3.24) 
OL/j m 

where 5j m is the Kronecker delta, then that 
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(hi)jj,j m (0,0) = (hi)jj >rn i(0,0) = (hi)j m: ji(0,0) = (h i)j m ,i n (p, 0) = 0, for all different 
j, m, I and n in 1^. Now, we set k = (tj m )j :Tn ^j z ^ in (I3.23D and we differentiate to get 

Oi)jm,Zn(0,0) = E tj 
a,/3,r),8 

which yields 

(/ii)jm,Zn(°>0) = T]i8j m 8i n + £i[6ji6 mn + S jn S m i], j, m, I, n G 1^, (3.25) 

where r/i and £j are real constants. 

We set in (pT2Tj) and (f3T22]) A ± = A 2 = C x = and B 1 = B 2 = B. Taking into account 
(pT20]) . (pT2l]l and (13T25D . we obtain, for {j,k} £ Si and {l,n} <£ S i} 

£?((Xi) jfc exptrS(X + r)<) = E (/^'^W 

= 0< E (£) iA,fc/ ^ 

= 9i(fi) jk , (3.26) 

and 

^((XOifcTOinexptrSCX + y)0 = E /f' ton<T [^^^ + ei(^V + ^ CT ^)] 

A,/i,!/,0"€/i-< 

= Vi(fi)jk,ln + £,i[(fi)jl,kn + (fi)jn,lkY (3.27) 

In the case where z is a separator, the system (|3. 17|) is replaced by 

£(exptr(,W^ + B X W^WU + dV^)) = ^(exptr^W^ + BiV^W£))£(exptrCiV^) 

= fi(Ax,B x )hi(Cx). 

and by the same reasoning, we also obtain equations (I3.26D and (I3.27p . but for j, k, I, n 
in h^. 

Now, to solve (|3.26p and (I3.27p . we introduce the functions 

ipi(B) = E(exptiBXi), ipi(B) = E(expirBYi). (3.28) 
Then (13.26P and (13.27P can be written as 

{^>i)jk^i = 0i(wPi)jk, (3.29) 

(<Pi)jk,ln^i = ViiVi^jkJn + €i[(<Pill>i)jl,kn + {Vi^jnjk] ' (3-30) 

Equation (I3.29P implies in particular that 

ipi = Mi) 6i . (3.31) 
Let <pi be such that, for B sufficiently close to zero, 

ifiipi = exp^j. (3.32) 
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Then (13.31H and (13. 32ft lead to the differential equation system defined for j, k, I, n £ Ih 
by 

Oi(<t>i)jk,ln + 0l{,4>i)jk{4>i)ln = £,i[(<t>i)jl,kn + (<t>i)jl{(f>i)kn + (4>i)jn,lk + (</>i) jn{<Pi)lk] 

+ riiK<f>i)jk,ln + {<fii)jk(<l>i)ln]- (3.33) 

The general solution of this system is 

(t>i{B) = PMB + D^+ai, 

where D £ V*,Xl £ ^> an d Oj G H are arbitrary constants and /3j is a constant depending 
on 9i, £j, r/j. This is in particular, justified by the fact that 

( B + D )fi = (-B + a) 

where A G V* . Now, if A = R*R £ P*, we denote A± opP = J R- 1 ( J R- 1 )* G P. Therefore, 
we use (|2.13p and a result in [2] due to Andersson and Wojnar which says that, for 6 = Z*Z, 
where Z G 7^ + , X = Z^e*^*)" 1 , to write 

4>i(B) = fit ln((B + D)^ opp ) M . + a- =/3i[lndet(S + D) opp ^ — lndet(-B + -D)^ opp J + a^, 

where a- is a constant. 
Thus 

= = H . and ^(B) 



(^Opp -jj.j (-L>^opp - 



I' i . 



We now observe that the fact that the t{j, i, j £ I are independent implies that the 
Xi, i £ I are independent. This is important for the calculation of the Laplace transform 
of X = ^^Xj. In fact, for each i £ I, by its very definition, Xi depends only on the 

iei 

such that i ^ j and on the i^- such that j, k G" and i ^ which are different from 
the ones on which any other component JQ/, i' ^ i depends. In other words, there exists a 
partition (rj)j e j of the set {t^, j, k £ 1} such that each Xi depends only on the tkj in Tj. 
Similarly the independence of the a^, i, j G / implies the independence of the Yi, i £ I. 
Finally denoting \ = {Pi^i, i £ 1} and x' = {PiO- ~ i G we obtain the Laplace 
transforms (/? of X = Xj and tpofY = Yi as 
iei iei 

Thus X ~ HW Xi(7 and Y" ~ HWv> where a = D x and a' = D x ' . Invoking Theorem 13.11 
i), we have necessarily D = e. This concludes the proof of the Theorem 13.21 □ 
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